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A numerical method is described for evaluating transverse spin correlations in the random phase 
approximation. Quantum spin-fluctuation corrections to sublattice magnetization are evaluated for 
the antiferromagnetic ground state of the half-filled Hubbard model in two and three dimensions 
in the whole U/t range. Extension to the case of defects in the AF is also discussed for spin 
vacancies and low— U impurities. In the U/t — * 00 limit, the vacancy-induced enhancement in the 
spin fluctuation correction is obtained for the spin-vacancy problem in two dimensions, for vacancy 
concentration up to the percolation threshold. For low-U impurities, the overall spin fluctuation 
correction is found to be strongly suppressed, although surprisingly spin fluctuations are locally 
enhanced at the low— U sites. 
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I. INTRODUCTION 

Transverse spin fluctuations are gapless, low-energy 
excitations in the broken-symmetry state of magnetic 
systems possessing continuous spin-rotational symme- 
try. Therefore at low temperatures they play an impor- 
tant role in diverse macroscopic properties such as exis- 
tence of long-range order, magnitude and temperature- 
dependence of the order parameter, Neel temperature, 
spin correlations etc. Specifically in the antiferromag- 
netic (AF) ground state of the half-filled Hubbard model 
transverse spin fluctuations are important both in two 
and three dimensions, where antiferromagnetic long- 
range order (AFLRO) exists at T = 0. In the strong- 
coupling limit (U/t — > 00), where spin fluctuations are 
strongest, they significantly reduce the zero-temperature 
AF order parameter in two dimensions to nearly 60% of 
the classical (HF) value. Similarly the Neel tem- 

perature in three dimensions is reduced to nearly 65% 
of the mean-field result T™ F = zJ/4 = 6t 2 /U, for the 
equivalent S — 1/2 quantum Heisenberg antiferromag- 
net (QHAF). @ 

Recently there has also been interest in spin fluctua- 
tions due to defects, disorder and vacancies in the quan- 
tum antiferromagnet. In the random-?/ model, where 
U is set to zero on a fraction / of sites, the lattice- 
averaged AF order parameter appears to be enhanced 
for small /, as seen in QMC calculations, ^ presumably 
due to an overall suppression of quantum spin fluctua- 
tions. On the other hand spin fluctuations are enhanced 
by strong disorder in the Hubbard model with random 
on-site energies. In the strong disorder regime, overlap 
of the two Hubbard bands leads to formation of essen- 
tially empty and doubly-occupied sites, which act like 
spin vacancies. [Q The problem of spin vacancies in the 
quantum antiferromagnet is also relevant to the electron- 
doped materials like Nd2_ x Ce x Cu04, where spin-dilution 
behavior is observed. [^],[)| While the problem of magnon 
energy renormalization due to spin vacancies has been 
addressed recently, JlO| [l5[] these methods are limited to 
the low-concentration limit, and the vacancy- induced en- 



hancement in transverse spin fluctuations has not been 
studied in the whole range of vacancy concentration. 

In this paper we describe a new method for eval- 
uating transverse spin correlations and quantum spin- 
fluctuation corrections about the HF-level broken- 
symmetry state, in terms of magnon mode energies and 
spectral functions obtained in the random phase ap- 
proximation (RPA). The method is applicable in the 
whole U /t range of interaction strength, and is illustrated 
with three applications involving the AF ground state of 
the half-filled Hubbard model — (i) the pure model in 
d = 2,3, (ii) spin vacancies in the strong coupling limit 
in d — 2, and (iii) \ow-U impurities in d = 2. This 
method for obtaining quantum correction to sublattice 
magnetization solely in terms of transverse spin correla- 
tions is parallel to the spin- wave-theory (SWT) approach, 
|]l6| , ^7| and differs from the method involving self-energy 
corrections. |l8| 

The RPA approach has been demonstrated earlier to 
properly interpolate between the weak and strong cou- 
pling limits for the spin-wave velocity. ]3]Jl8|| By go- 
ing beyond the RPA level within a systematic inverse- 
degeneracy expansion scheme, which preserves the spin- 
rotational symmetry and hence the Goldstone mode or- 
der by order, it was also shown that in the strong coupling 
limit identical results are obtained for all quantum cor- 
rections, order by order, as from the SWT approach for 
the QHAF. Q A renormalized RPA approach has also 
been devised recently to obtain the magnetic phase di- 
agram for the three dimensional Hubbard model in the 
whole U/t range, and the Tn vs. U behaviour was shown 
to properly interpolate between the weak and strong cou- 
pling limits. Q 



II. TRANSVERSE SPIN CORRELATIONS 



The method is based on a convenient way to perform 
the frequency integral in order to obtain spin correla- 
tions from spin propagators, and we illustrate it here 
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for transverse spin correlations. We write the time- 
ordered, transverse spin propagator for sites i and j, 
{y G \T[Si(t)Sf(t')]\V G ) at the RPA level in frequency 
space as, 
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where \<p n (uj)} and X n (u>) are the eigenvectors and 
eigenvalues of the matrix. Here [x°(ui)]ij — 

iJ(du'/2Tr)G] j (u')G^(cj' - w) is the zeroth-order, 
antiparallel-spin particle-hole propagator, evaluated in 
the self-consistent, broken-symmetry state from the HF 
Green's functions G <T (u>). Spin correlations are then ob- 
tained from, 
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where the collective mode energies u n are obtained from 
1 — UX n (u> n ) = 0, and X n (u>) has been Taylor-expanded 
as X n (u>) ~ X n (Lu n ) + (u> — Lo n )(dX n / 1 doj) Un near the mode 
energies to obtain the residues. For convergence, the re- 
tarded (advanced) part of the time-ordered propagator 
X h , having pole below (above) the real-w axis, is to be 
taken for t' < t (t' > t). The frequency integral is con- 
veniently replaced by an appropriate contour integral in 
the lower or upper half-plane in the complex-w space for 
these two cases, respectively, which results in Eq. (2). 



III. HUBBARD MODEL IN D = 2,3 

We first illustrate this method for the half-filled Hub- 
bard model in two and three dimensions on square and 
simple-cubic lattices, respectively. In this case it is con- 
venient to use the two-sublattice representation due to 
translational symmetry and we work with the 2x2 ma- 
trix [x m momentum space, which is given in terms 
of eigensolutions of the HF Hamiltonian matrix. The 
fc-summation is performed numerically using a momen- 
tum grid with Ak = 0.1 and 0.05, in three and two di- 
mensions, respectively. 

Equal-time, same-site transverse spin correlations are 
then obtained from Eq. (2) by summing over the dif- 
ferent q modes, using a momentum grid with Aq = 0.3 
and 0.1 in three and two dimensions, respectively. We 
consider t' — > t~, so that the retarded part is used, with 
positive mode energies. From spin-sublattice symmetry, 
correlations on A and B sublattice sites are related via 
(S + S~)a = (S~S+)b- Thus the transverse spin correla- 
tions are obtained from magnon amplitudes on A and B 
sublattices, and from Eq. (2) we have 
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From the commutation relation [S ,+ ,5 f_ ] = 2S Z , the 
difference (S + S~ — S~S + }rpa, of transverse spin corre- 
lations evaluated at the RPA level, should be identically 
equal to (2S' z )hf- This is becasue both the RPA and HF 
approximations are 0(1) within the inverse-degeneracy 
expansion scheme H in powers of 1/M {M is the num- 
ber of orbitals per site), and therefore become exact in 
the limit N — * oo, when all corrections of order 1/Af 
or higher vanish. This is indeed confirmed as shown in 
Figs. 1 and 2. The deviation at small U is because of 
the neglect in Eq. (2) of the contribution from the single 
particle excitations across the charge gap, arising from 
the imaginary part of ~x°{oj) in Eq. (1). 

The sum (S + S~ + S~S + )rpa yields a measure of 
transverse spin fluctuations about the HF state, and in 
the strong coupling for spin S, one obtains {S + S~ + 

S-S+) RPA = (2S) J2 q V^l ~7 9 2 - IB Using the iden- 
tity, (S Z S Z ) = S(S+l)-(S+S- + S-S+)/2 in this limit, 
the sublattice magnetization m = {2S Z ) is then obtained 
from, 



(S z ) = S 



_ 1 / (s+s- + S~S+) _ 
S\ 2S 



1/2 



(4) 



To order 1/2S, this yields the correction to the sublat- 
tice magnetizationof (2S , )- 1 ^ <3 [(l- 7 2)- 1 /2_ 1 ] = .156 

and 0.393 for S = 1/2 in three and two dimensions, re- 
spectively This same result at one loop level was also ob- 
tained from a different approach in terms of the electronic 
spectral weight transfer, [p|j4[| and is in exact agreement 
with the SWT result. @|j 

As (S' z )hf is the maximum (classical) spin polarization 
in the z-direction, and therefore also the maximum eigen- 
value of the local S z operator, therefore for arbitrary U, 
the HF magnitude (S z )hf plays the role of the effective 
spin quantum number, S. The sublattice magnetization 
m is therefore obtained from m = tohf - Stusf, where the 
first-order, quantum spin-fluctuation correction Stusf is 
obtained from Eq. (4) with S — (S z )uf, 
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In the strong coupling limit U/t 



oo, (S+S- - 

S-S+) RP a = (2S Z ) HF = 1 for S = 1/2, so that 
the spin-fluctuation correction simplifies to, Smg-p = 
2(S~ S+) R pa- For a site on the B sublattice, where 
(S+S- - S^S^rpa = (2S z )hf = -1, we have 5m SF = 
2(S+S-) RPA . 
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FIG. 1. The sublattice magnetization m vs. U in 
two dimensions (diamonds), along with the HF results 
from (i) the self-consistency condition (dashed), and (ii) 
<25 z )hp = (S+S- - S-S+}rpa (plus). 
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FIG. 3. The sublattice magnetization m vs. U from our 
RPA spin-fluctuation approach (solid line) compared with 
earlier results — VMC [19] (diamonds), SE [18] (dash-dot), 
QMC [20] (errorbars), and GLSWA [22] (dash). Also shown 
is the HF result (dotted). The arrows denote the asymptotic 
results of (a) [2] and (b) [1]. 
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FIG. 2. The sublattice magnetization m vs. U in 
three dimensions (diamonds), along with the HF results 
from (i) the self-consistency condition (dashed), and (ii) 
<2S z )hf = (S+S- - S-S+Upa (plus). 

The [/-dependence of sublattice magnetization m is 
shown in Figs. 1, and 2 for d = 2 and d — 3, respec- 
tively. In both cases it interpolates properly between the 
weak and strong coupling limits, approaching the SWT 
results 0.607 and 0.844, respectively, as U/t — > oo. A 
comparison of the m vs. U behaviour with earlier re- 
sults is presented in Fig. 3 for the well studied d = 2 
case. Earlier studies have employed a variety of meth- 
ods including the variational Monte Carlo (VMC), JToj ] 
self-energy corrections (SE), |18| quantum Monte Carlo 
(QMC), [^0| functional- integral schemes, |2lJ] the general- 
ized linear spin- wave approximation (GLSWA), p2j and 
mapping of low-energy excitations to those of a QHAF 
with [/-dependent, extended-range spin couplings. p3[ 



In addition to the two-sublattice basis, we have also 
used the full site representation in the strong coupling 
limit, in order to illustrate the scaling of the quantum 
correction with system size. Here the x° matrix is evalu- 
ated and diagonalized in the site basis for finite lattices. 
Results for lattice sizes with 8 < L < 16 are shown in Fig. 
4. A quadratic least-square fit is used to extrapolate to 
infinite system size, which yields Sma,-p(l/L — > 0) = 0.39, 
in agreement with the result from Eq. (4). 



IV. SPIN VACANCIES 

The site representation has also been used to obtain 
transverse spin fluctuations for the problem of spin va- 
cancies in the AF in the limit U/t — > oo. As mentioned 
already this method is applicable in the whole range of 
vacancy concentration, and allows determination of the 
critical vacancy concentration at which the AF order 
parameter vanishes. For the vacancy problem we con- 
sider the following Hamiltonian on a square lattice with 
nearest-neighbor (NN) hopping, 
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where the hopping terms tij = if sites i or j are vacancy 
sites, and £y = t otherwise. Thus, for a vacancy on site i, 
all hopping terms t^ connecting i to its NN sites j are set 
to zero. The vacancy site is thus completely decoupled 
from the system. Half-filling is retained by having one 
fermion per remaining site. We consider the U/t — ► oo 
limit, where the local moments are fully saturated, and 
the vacancy problem becomes identical to the 
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FIG. 4. The spin-fluctuation correction to sublattice mag- 
netization SmsF vs. 1/L, for L — 8, 10, 12, 14, 16. 

spin-vacancy problem in the QHAF. This is also equiv- 
alent to the problem of non-magnetic impurities in the 
AF in the limit of the impurity potential V — * oo. [ p^Jl5| 
The structure of the X (w) matrix in the host AF, and 
the modification introduced by spin vacancies has been 
considered earlier in the context of static impurities. Jjl| 
As the Goldstone mode is preserved, it is convenient to 
work with the matrix K = U[l — U\ a ], inverse of which 
yields the spin propagator in Eq. (1) near the mode en- 
ergies. When expressed in units of U 2 t 2 /A 3 — 2 J, where 
2A = mnpU is the AF gap parameter, and J — 4t 2 /U, 
it has the following simple structure for the host AF: 

Kii = l+o; i in A sublattice 
= l — o; i in B sublattice 

1/4 jisNNofz. (7) 



K°. 



Vacancies introduce a perturbation in K due to ab- 
sence of hopping between the vacancy and NN sites, and 
we take 5K = —U 2 Sx° to refer to this vacancy-induced 
perturbation, so that K = K a + 5K. If site i is a vacancy 
site and j the NN sites, then the matrix elements of 5K 
are, 



SKa = 6K M = SKjj = -1/4 , 



(8) 



and the magnitude of Ka is irrelevant since the vacancy 
site i is decoupled. Thus Kij = 0, reflecting the decou- 
pling of the vacancy, and the static part of the diagonal 
matrix elements Kjj on NN sites are reduced by 1/4. For 
n vacancies on NN sites, the static part is 1 — n/A. This 
ensures that the Goldstone mode is preserved. Thus if 
a spin on site j were surrounded by a maximum of four 
vacancies on NN sites, then the static part vanishes, and 
Kjj = ±lu, representing an isolated spin, which yields a 
1/lu pole in the transverse spin propagator. This is an 
isolated single-spin cluster, and with increasing vacancy 
concentration, larger isolated spin clusters are formed. 
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FIG. 5. The spin-fluctuation correction to sublattice mag- 
netization SmsF vs. x, for L = 10 (diamonds), 12 (plus), and 
14 (squares). 

As these are decoupled from the remaining system, their 
spin-fluctuation contributions are not included, as the 
OP vanishes for finite spin clusters. When x exceeds 
the percolation threshold ~ 0.4, the fraction of macro- 
scopically large spin clusters in the system vanishes, and 
therefore no AFLRO is possible for x > 0.4. 

For a given vacancy concentration and system size, 
the appropriate number of vacancies are placed randomly 
across the lattice, and the matrix K constructed accord- 
ingly. Exact diagonalization of K is carried out, and the 
eigensolutions are used to compute the transverse spin 
correlations from Eq. (2). The quantum, spin-fluctuation 
correction is then obtained from the strong-coupling limit 
of Eq. (5). The transverse spin correlation (S~S + ) is 
averaged over all spins within the A sublattice. As men- 
tioned already, only the spins in the macroscopic cluster 
spanning the whole lattice are considered, and contribu- 
tions from spins in isolated spin clusters are excluded. 
Configuration averaging over several realizations of the 
vacancy distribution is also carried out. 

The quantum spin-fluctuation correction vs. vacancy 
concentration is shown in Fig. 5 for three lattice sizes, 
L = 10, 12, 14. Best fits are obtained with an expres- 
sion including a cubic term, SmsF — a + fix + jx 3 . For 
the three lattice sizes a =0.236, 0.260 and 0.278, respec- 
tively, and as shown in Fig. 3, it extrapolates to 0.39 as 
1/L — » 0. The coefficient of the linear term is found to 
be nearly independent of system size, (3 « 0.42. And the 
cubic term 7 takes values approximately 3.6, 4.0, and 4.4 
for the three lattice sizes, and extrapolates to about 6.5 
as 1/L — > 0. With these coefficients, the spin-fluctuation 
correction Stusf is nearly 1 for x = 0.4. Beyond x = 0.4, 
the percolation limit, there is no single, macroscopically 
large spin cluster left in the system. Therefore the point 
where the AF order parameter vanishes and AFLRO is 
destroyed due to spin fluctuations nearly coincides with 
the percolation threshold. This is in agreement with 
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results from series expansion p40 and quantum Monte 
Carlo simulations Pa] of the QHAF with spin vacancies. 



V. RANDOM-?/ MODEL 

We now consider a quenched impurity model with a 
random distribution of impurity sites characterized by 
a local Coulomb interaction U' ^ U for the host sites. 
With H and I referring to the sets of host and impu- 
rity sites respectively, we consider the following Hubbard 
model in the particle-hole symmetric form at half-filling 
and on a square lattice, 

H = -t J2 ( a l a j* + h - c -) + u IDfat ~ \ ){n ^ - b 
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The motivations for studying this impurity model are 
threefold. In view of the observed enhancement of mag- 
netic order at low concentration of impurities, || we shall 
analyze the suppression of quantum spin fluctuations to 
examine whether this is due to a local suppression at the 
low-C/ sites. The RPA evaluation of transverse spin cor- 
relations is also extended to the case of site-dependent 
interactions. Furthermore, at half-filling this model also 
provides a simplistic representation for magnetic impu- 
rity doping in an AF. This may appear contradictory in 
view of the apparently nonmagnetic (U' « 0) nature of 
the impurity sites. However, this feature is expressed 
only away from half-filling. || 

The atomic limit t = provides a convenient starting 
point for further discussions. In the particle- hole sym- 
metric form of Eq. (9), since not only local interaction 
terms, but the on-site energy terms are also modified 
(from —U/2 to —U'/2) at the impurity sites, therefore the 
energy levels for added hole and particle are (—U/2, U/2) 
and (—U'/2,U'/2) for host and impurity sites respec- 
tively. To order t 2 , this impurity model therefore canon- 
ically maps to the following S — 1/2 Heisenberg model 



H c tf — J 
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where in the first term J = At 2 /U is the conventional 
exchange coupling between neighboring host spins, and 
J' = 8t 2 /(U + U') is the exchange coupling between im- 
purity spins and neighboring host spins. In writing Eq. 
(10) we have assumed the dilute impurity limit, and dis- 
counted the possibility of two impurity spins occupying 
NN positions, in which case the impurity- impurity ex- 
change coupling will be At 2 /U 1 . 

Therefore, in the strong correlation limit, this random- 
ly model also describes magnetic impurities in the AF 



within an impurity-spin model. The magnetic-impurity 
doping is characterized by identical impurity and host 
spins (S 1 = S) , but with different impurity- host exchange 
coupling (J' / J). Both cases J' > J or J' < J are possi- 
ble, and in this paper we have considered the two cases: 
(i) U' « U so that J' « 2 J, and U' = 3U so that 
J' = J/2. While this model is easily generalized to other 
magnetic-impurity models represented by locally modi- 
fied impurity- host hopping terms t' ^ t, and/or different 
impurity energy levels, in fact, the essential features are 
already contained here as the impurity exchange coupling 
J' is the relevant quantity in determining the spin fluc- 
tuation behavior. 



RPA with site-dependent interaction 

We recast the RPA expression for the transverse spin 
propagator in a form suitable for site-dependent inter- 
actions. In terms of a diagonal interaction matrix [U], 
with elements [U}u — Ui, the local Coulomb interaction 
at site i, the time-ordered transverse spin propagator at 
the RPA level can be rewritten, after simple matrix ma- 
nipulations, as 

[Y -+y-JVL- 1 1 fin 

where [•A(w)] = [U] — [U][x°{u> )][U] is a symmetric ma- 
trix. As [U] is non-singular, the singularities in [\ h (w)], 
which yield the magnon modes, are then given com- 
pletely by the vanishing of the eigenvalues of the ma- 
trix [A]. In terms of A n and \(j> n ), the eigenvalues and 
eigenvectors of the matrix [A], we have [A(w)] _1 = 
A„(ti>) _1 |(/)„(a;)}(|(/> n (a;)|, so that the magnon- mode 
energies aj n are then given by A ra (w ra ) = 0, and in anal- 
ogy with Eq. (2) the transverse spin correlations are 
obtained from, 



(s-(t)s+(t')U PA = ±J2 



<l>n(UnWn{Un) 

{d\ n /dw) Un 
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As we are interested in the dilute behaviour, we exam- 
ine the correction to sublattice magnetization due to two 
impurities, one on each sublattice for symmetry. Since 
m = m^F — SmsF, corrections to both mnF and SmsF are 
expressed in the dilute limit (impurity concentration x) 
as tohf = "t-hf — ctiiFX, and <5msF = <5mg F — asF^c. The 
overall m(x) behavior therefore depends on the relative 
magnitude of the coefficients anF and ctsF- For U = 10 
and U' = 2 we find, at the HF level, that m^ F = 0.93 
and ohf = 0.19, both quite independent of system size. 

The (site-averaged) spin fluctuation correction t>rnsF 
is obtained from Eq. (5) with and without impurities, 
and the impurity contribution extracted. For U' = 2 we 
find a net reduction in Suisf- Divided by 2/N, the impu- 
rity concentration, this yields the coefficient asF defined 



5 



above, and also the per-impurity contribution to the to- 
tal spin fluctuation correction over the whole lattice. The 
spin-fluctuation correction <5msF for the pure case, and 
the per-impurity contribution (o!sf) are shown in Fig. 6 
for different lattice sizes, along with least-square fits. It 
is seen that in the infinite size limit, the per-impurity 
reduction is nearly 0.2, which is more than half of the 
correction per site in the pure case (0.35). Thus, there is 
a substantial reduction in the averaged spin fluctuation 
correction due to the low— U impurities. 

As the two coefficients «hf an d «sf are very nearly 
the same, the sublattice magnetization 

m{x) = m(0) — (q!hf — osf)x (13) 

shows negligible concentration dependence. Thus the 
(relatively small) reduction in the HF value due to the 
low-iJ impurities is almost fully compensated by the (rel- 
atively substantial) reduction in the spin fluctuation cor- 
rection. To first order in x, we thus find that there is 
nearly no loss of AF order due to the law-U impurities. 
As mentioned earlier, even a slight enhancement in the 
AF order was recently seen for the case U' = 0. || 

We next examine the site-dependence of the local spin 
fluctuation corrections Sm l SF near the impurities. Table 
I shows that spin fluctuation is actually enhanced on the 
low-U impurity sites. The suppression of 5m l SF in the 
vicinity more than compensates for this local enhance- 
ment, resulting in an overall reduction on the average. 
On the other hand, for U' = 30, we find that the correc- 
tion is suppressed on the high-L 7 impurity site, while it 
is enhanced on the average. Thus, to summarize, when 
the impurity spin is coupled more strongly (weakly), the 
spin-fluctuation correction is enhanced (suppressed) lo- 
cally at the impurity site, but the average correction to 
sublattice magnetization is suppressed (enhanced). 
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FIG. 6. The per-impurity reduction in the spin-fluctuation 
correction to sublattice magnetization SmsF vs. l/L (plus), 
for finite lattices with L = 8, 10, 12, 14, 16. Also shown is the 
correction per site for the pure system (diamonds), indicating 
significant relative impurity contribution. 

This local enhancement can be understood in terms 
of the correlations (S^S^~) as follows. Since the im- 
purity spin is more strongly coupled to the neighboring 
spins, the NN matrix elements are enhanced. This 
puts more magnon amplitude </> on the impurity site, 
so that from Eq. (2) the transverse spin correlations 
{S~S^~), and therefore the spin-fluctuation correction, 
are enhanced for low-U impurities. The overall decrease 
in the averaged fluctuation correction, however, is due to 
the stiffening of the magnon spectrum in the important 
low-energy sector, following from the increased average 
spin coupling. 



TABLE I. The local spin-fluctuation corrections Sml t for a 16 x 16 system (U = 10), with two impurities at (11,4) and 
(4,14). For the two cases U' = 2 (J' ~ 2J) and U' = 30 (J' = J/2), quantum corrections are enhanced/suppressed locally at 
the impurity sites (indicated in boldfaces), but are suppressed/enhanced on the average. 
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.254 
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.253 


.265 
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.265 
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.254 


.253 
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.253 


.261 


.248 


.158 


.248 


.261 


.253 


.253 


15 


.254 


.253 


.264 


.248 


.265 


.253 


.253 


.253 


16 


.253 


.253 


.253 


.260 


.253 


.253 


.253 


.253 



G 



While these results also follow from the impurity-spin 
picture, the small charge gap at the impurity site does 
have an impact on the magnon spectrum. Within the 
localized spin picture, the highest-energy magnon mode 
corresponding to a local spin deviation at the impurity 
site would cost energy 4 x J'/2 = 16t 2 /(U + U'). How- 
ever, the highest energy in the magnon spectrum is actu- 
ally seen to be 1.05i, which is substantially smaller than 
2,1' = 1.33i. This shows the compression effect of the low 
charge gap (2.75i) on the magnon spectrum. |2q| 

In conclusion, using a convenient numerical method 
for evaluating transverse spin correlations at the RPA 
level, quantum spin-fluctuation corrections to sublattice 
magnetization are obtained for the half-filled Hubbard 
model in the whole U/t range. Results in two and three 
dimensions are shown to interpolate properly between 
both the weak and strong correlation limits, and ap- 
proach the SWT results as U/t — > oo. The method is 
readily extended to other situations of interest involving 
defects in the AF, such as vacancies/impurities/disorder. 
Numerical diagonalization for finite lattices, along with 
finite-size scaling tested with the pure Hubbard model, 
allows for exact treatment of defects at the RPA level. 
This is illustrated with a study of the defect-induced 
enhancement / suppression in transverse spin fluctuations 
for spin vacancies and low-J7 impurities in two dimen- 
sions. While the quantum spin fluctuation correction to 
sublattice magnetization is sharply enhanced by spin va- 
cancies, it is strongly suppressed by the low-f impuri- 
ties, although the fluctuation correction is enhanced at 
the low— U sites. 
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